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We have recently argued that if one introduces a relational time in quantum mechanics and 
quantum gravity, the resulting quantum theory is such that pure states evolve into mixed states. 
The rate at which states decohere depends on the energy of the states. There is therefore the 
question of how this can be reconciled with Galilean invariance. More generally, since the relational 
description is based on objects that are not Dirac observables, the issue of covariance is of importance 
in the formalism as a whole. In this note we work out an explicit example of a totally constrained, 
generally covariant system of non-relativistic particles that shows that the formula for the relational 
conditional probability is a Galilean scalar and therefore the decoherence rate is invariant. 

I. INTRODUCTION 

We have recently introduced a new technique for discretizing physical theories P]. When applied to general rela- 
tivity it yields a discrete theory that is constraint-free yet it approximates well the continuum theory under certain 
circumstances 0, |^ ■ The lack of constraints allows to tackle some of the fundamental open problems of canonical 
quantum gravity. For instance one can introduce a relational time a la Page-Wootters '?]. That is, one pro- 

motes all quantities in the theory to quantum operators and chooses one that is called "clock" and then one computes 
conditional probabilities for the other variables to take given values when the "clock" variable shows a certain "time" . 
The resulting quantum theory approximates ordinary quantum mechanics well when the clock variable chosen behaves 
in a semi-classical fashion with small quantum fluctuations. If one chooses as clock a variable that is in a quantum 
regime, the resulting theory is still valid but it will not resemble ordinary quantum mechanics. 

We have also argued that, due to the fact that one cannot have a perfectly classical clock in nature, the resulting 
theory will have small but non- vanishing departures from ordinary quantum mechanics. In particular a pure state 
does not remain pure forever but evolves into a mixed state. 

Since one is approximating a constrained continuum theory with a discrete theory that is unconstrained, the 
resulting relational discrete theory is formulated in terms of variables that are observables for the discrete theory. But 
they are not necessarily the discrete counterpart of Dirac observables of the continuum theory. Therefore the issue 
of how to reconcile the predictions of the discrete relational theory with the covariance of the continuum theory is of 
importance. In particular, the conditional probabilities must remain invariant when one changes coordinates and both 
the clock variable and the observed variable change values. To tackle the covariance problem in complete generality 
is beyond the scope of this paper. What we intend to do here is to analyze a simple model where calculations can be 
worked out concretely and in particular to probe the following issue. Since the prediction for the time of decoherence 
of pure states results in a formula that involves the energy of the states, it is may not be immediately apparent in 
what sense it is Galilean invariant. We would like to discuss in a simple model how to interpret the formula in a way 
that the invariance is manifest. 

The organization of this paper is as follows. In section II we present the model we will study and in section III we 
will show the emergence of Galilean invariance. We end with a discussion. 

II. THE MODEL 

We consider the following model. It consists of two non-interacting particles moving in separate potentials in 1 + 1 
dimensions. One of the particles we will assume is much more massive than the other and it will determine the 
variable we choose as a clock. The other particle we will call the "system" particle. The potential affecting the clock 
particle will be a constant force field. We will assume the particle is far away from the turning point, since in this 
regime we know our discrete approach approximates well the continuum I3| . For the system particle we will assume it 
behaves quantum mechanically and is in a potential that gives rise to bound states. As is well known 8], the best way 
to understand Galilean transformations in quantum mechanics is to study them as a limit of Lorentz transformations. 
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We will therefore choose Lorentz invariant and reparameterization invariant action for the particles, 
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where g°,(7 are the space-time coordinates of the "clock" particle and '/>°,(/' are the space-time coordinates of the 
"system" particle. We have kept the speed of light explicit in order to consider later the non-relativistic limit. We 
start by assuming that the time-like coordinates of both particles can be synchronized (since we will work in the 
Newtonian limit this poses no conceptual problem) = 0°, and the reference system has been chosen such that 
(7° (7 and (jfi » and also Mc^ U{q^,q) and mc^ 3> V{(j)°,(j)) (non-relativistic limit). For concreteness we 
assume that in the reference frame given, U{q^, q) = aq and V{(l)°,4>) = V{(j)) and the latter has bound states. With 
these assumptions the action becomes. 
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where A is a Lagrange multiplier associated with the constraint that imposes the synchronization. It is immediate 
to see that if one chooses q^ = (jp = ct one will obtain the ordinary action for two non-relativistic particles with t 
the ordinary non-relativistic time. We will not do this here, since we are interested in handling a totally constrained 
system, since it is in such systems where the introduction of a relational time is meaningful since they have no preferred 
notion of time. 

To understand better the constraint structure of the theory, we will rewrite the action in first-order form. We define 

the canonical momenta. 
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From where we can get two constraints, in addition to the one we had before (jp — q^ = 0, 
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If we rearrange the latter two constraints into their sum and difference, 
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one readily sees that the last constraint together with q^ — (fp = are second class, whereas they both commute with 
(PJ. One imposes the second class constraints strongly and is left with a theory with one constraint, whose action is. 
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where we introduced the shorthand po = pq + ttq. This action is very natural for the system under study (in fact we 
could have started the calculation simply by considering this action from the outset). 
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III. GALILEAN INVARIANCE 



To probe the invariance of the decoherence effect of interest, we would hke to study two different situations. One in 
which the system particle is a potential V{(l)) and another in which the potential is of the form V{4> — (3q^), this will 
represent a system that is bound by a potential around some minimum that is fixed or that is moving, respectively. 
This corresponds to adopting the active point of view of the Galilean transformation. We will do this with our 
consistent discretization techniques We refer the reader to our previous papers for details on the technique. 

We start by discretizing the action in the first of the two cases of interest. The integral in the action becomes replaced 
by a discrete sum S = L{n, n + 1) and we absorb the time interval e = t„+i — r„ and where. 



n+1) = Pan {cin+l - 9°) + Pn (<7n+l - qn) + T^n (</'n+l " " 
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We now implement the canonical transformation that materializes the time evolution between instant n and n+1 
with the Lagrangian —L(n,ji + 1) playing the role of generating function of a type I canonical transformation, 
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The system has constraints and we will use them to eliminate some of the variables and yield a system of evolution 
equations in a more explicit form. The resulting system is the following, 

(28) 
(29) 

(30) 
(31) 



Hn+1 


= 9^ 




pq" 

n+l 


n 




Qn+l 


= Qn 4 


pi 

M 


pq 

n+l 


= pi - 

n 





4 



„+i = 0„ + 7V„:^, (32) 
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The last equation determines the Lagrange muhiphcr Nn- To see this, we first rewrite it entirely in terms of 
variables at n + 1 , 

g^^-iV^ (i^ ^ n0»+i-jV„P„^^,) ^ (P„\,)^ 

Since we are ultimately interested in studying the system in a regime close to the continuum limit, we make the 
assumption that the lapse Nn is small and expand the term involving the potential to first order in Nn, 
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We can now solve explicitly for the Lagrange multiplier, 

^n+l , ^n+lTTi 



^« = + ^^'(-^rO Cn+1 (37) 

\ Aic TOc ' 
Where C„+i is the constraint of the continuum theory discretized as if all variables were at n + 1, 

uq" > AT i 1"+^ I (-^n+l)^ , , ^i't'n+l) , (^n+l)^ ,gnx 

C„+i = P„\, + A/c + «^ + + mc + + . (38) 

We now assume that a ^ V {(/)). This is due to the fact that we are assuming the clock to be classical and large and 
a is therefore associated with a macroscopic force whereas V{(f>) is the potential in which the system is bound, and 
the latter is microscopic in nature. With this assumption we make sure there are no singularities in the computation 
of the Lagrange multiplier. Recall that the discrete description departs from the continuum one close to the turning 
point of the orbit. 

One can now substitute the expression of the Lagrange multiplier in the evolution equations. The resulting system 
of equations can be viewed as a canonical transformation between instant n and instant n + 1 for the remaining 
variables of the problem. The next step consists in quantizing the system by representing the discrete evolution 
through a unitary operator, i.e. £^ = Wzl^^j^U where the z*'s are all the phase space variables of the problem. All 
these calculations can be worked out explicitly for a simple system like the one we are considering, we will not show 
all the details here for reasons of space, the reader can see similar treatments in 

Since we are interested in the continuum limit, a shortcut can be taken by considering the Hamiltonian associated 
with the unitary transformation U = e*^ . The Hamiltonian is obtained by taking the logarithm of the unitary 
operator as a power series. This power series is convergent at all points in phase space except for a small region 
around the turning point of the orbit of the clock system. The Hamiltonian is obviously conserved upon evolution 
(except at the turning point). The first term in the expansion of the Hamiltonian is. 
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and to simplify notation, from now on we call — p„ and P^^ = pon- 

For the quantization we consider wavefunctions il}n{<f ,<It4>) forming a Hilbert space at the "instant" n. Isomor- 
phic Hilbert spaces exist at all other discrete instants. With the Hamiltonian we will study the evolution operator 
U{n,nfj) — Q^H{n-no) g^^^ j^g action on the states, ■0„((7°, g, (/)) — lJ{n,no)ipno{'f t1^4>)^ and we are working in the 
Schrodinger representation. The explicit form of the quantum Hamiltonian is. 



Mc / - - \ 2 
H = —r(po + Hi+H2) . (40) 
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It is to be noted that the expression in parenthesis is the constraint that one has in the continuum theory. In the 
consistent discretization approach the constraint of the continuum theory is not enforced exactly (what is enforced 
is equation H34|l which corresponds to the constraint of the continuum theory but with the momenta evaluated one 
instant after the configuration variables). In the continuum limit, it nevertheless is enforced quite approximately and 
therefore the norm of H is going to be small. 

We consider a quantum state in which the clock has a semiclassical behavior, so it is described by a coherent state 
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where cti is the dispersion in the variable q and ctq is the dispersion in the variable We have also assumed that 
E ^ \po + E\ < H2 >no- The first inequality is in order to be in the continuum limit. The second inequality is in 
order to simplify calculations, and implies that we are accepting as "continuum limit" a regime where the constraint 
of the continuum theory is well enforced with respect to the scale of energies of relevance for the "clock" system, but 
the error in enforcement is large with respect to the energies of the system under study. It would be desirable to 
extend the results of this paper to regimes that approximate even further the continuum theory, but the calculations 
would be more involved. 

The fundamental equation to be studied is the conditional probability. 
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with P(f,^qO is the projector onto the eigenstate labeled by the values 0, q'^ and p^o, pq, the density matrices associated 
with the state '^no- From now on we will use natural units where c = h = 1. 

Let us analyze the denominator of this expression. Taking the trace on the 0, q spaces by integrating, we get. 

Den = ^Tr (u^n)PqoU{n)pqO x pg x pA (44) 
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and the term involving from U cancels with that of W since commutes with P^o . We have also replaced Hi 
by E and p by p since the trace implies taking the expectation value of quantities depending on q and p. Since PqO 
represents a state very peaked at < po > — Po and < if >— and since |po + ^ < H2 >, we have that. 
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where Pn[q^) = Pqapn,q° = PqoU{n)pqoU'' (n) represents the wavepacket of a "free particle" which evolves with the 
effective Hamiltonian 



ap 

u 



(46) 



It is instructive to realize that one can write. 
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where 
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which shows that the clock "displays a time" in the neighborhood of (f when we are at the level n of the discrete 
theory. We have defined tmax('T-), the most likely value of the clock "time" for a given n level in the discrete theory, 
and we have chosen the clock in such a way that tmax grows linearly with n. 
The width of the packet grows with n as, 
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We now should introduce some relevant scales. We will assume the characteristic mass of the clock system is about 
a kilogram. The potential of the clock system is characterized by the macroscopic constant a, which we will assume 



is of the order of 10 Newton, which in natural units corresponds to a 
(To ~ 10^^"s ^ Im, then. 
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and we have assumed n/M ^ 10~^ so we are in a non-relativistic regime. 

As we discussed in the sums that appear in the numerator and denominator for the conditional probability 
should be large enough to involve the complete evolution of interest for the system, but they should not be infinite, 
otherwise one gets an indeterminate quotient of two diverging quantities for the conditional probability. Given the 
value computed above for the quantity multiplying (n — no)^, it is natural to bound the value of n — rip ^ 10^^, that 
is we assume that the sums go from no to a maximum value N ^ 10^^, let's say TV ^ 10^*^, otherwise the packet 
representing the clock will spread too much and we would be out of the semiclassical regime. Notice that we also 
have that E > lO^^m"^, and recalling that \po + E\ has to be smaller than E, and choosing it to be lO^^m"^ yields 



10-*' 



s ^ 3 hours, which is a reasonable number. Summarizing, by bounding the number of steps we find that the 
denominator is a quantity of order unity. Its precise value is not of great interest, since we can choose it by fixing the 
normalization of the probability. 
Let us analyze the numerator, 
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and observing that the projector is independent of q, and one can therefore substitute Hi by its expectation value E. 
Using now that the clock is semiclassical and po -'r E ^ < H2 > to neglect terms quadratic in H2 we can write, 
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where we have replaced po with po since the clock is approximately classical and its energy dominates in po, and we 
have separated the expression into two pieces, one dependent on the variable and one dependent on the qo variable 
to make more explicit the separation between clock and system. 

The last trace divided by Den can be written as VnisP) and satisfies that 'Ylm '^nifp) = 1- 
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Following the discussion in |y] , in order to make contact with ordinary quantum mechanics we assume the spacing 
in n is small compared with the values of n and introduce a continuous variable v = nt. We choose e such that 



ap 



(54) 



so we have that e < Im with the choice of scales we made for the problem. We choose no = and we can then write 
a good continuum limit approximation for 'P,i{q°), as in 0, 
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I = V, so we can write, 



, 1 (q"? 



4 



4r74 4(po + EY 



(55) 



(56) 



and with imax('^) = f , 
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and this density matrix satisfies a Schrodinger equation modified due to the fact that we are considering a quantum 
clock as shown in detail in |g, 



dp2 



t[H2,P2]~<T{q'')[H2,[H2,P2]], 



(59) 



with <T{q'^) — daQ[q^) / dq^ . This expression is just the first two terms in a power series in terms of the dispersion of 
the quantum clock, which for realistic systems is a very small quantity. 

To try to get a handle on a rough value for this quantity in the case of a realistic system, we note that the 
macroscopic clock particle is subject to decoherence due to interaction with the environment. If we characterize such 
decoherence by a time t_D, we have, 
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If ~ Is 10^°m, which is a rather large decoherence time for a macroscopic system, then a ~ In 
reference |lOj we have estimated theoretical limits as to how small a dispersion is attainable with optimal realistic 
clocks. 

In order to study the Galilean covariance of the conditional probability, the procedure is simple. We have to repeat 
the calculation assuming a boost with velocity —j3 has been performed on the system 2 respect to the system 1, in 
such a way that the potential it now sees is of the form V{(j) — (iq^)- For instance, the system 2 can be an electron in 
a central potential given by a nucleus. The relational analysis goes along exactly as before, with two differences. The 
Hamiltonian for the second system becomes. 
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and the initial state of the system is given by p'^o ^ p'^^ Uq {pqo x p^) Lf^ with 
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In other words, the initial state is the one corresponding to the Galilean boost Uq to the original state Notice 
however that in traditional treatments of Galilean invariance in quantum mechanics the variable that we here take 
as (f is a classical parameter t. Our treatment can be considered a relational generalization of the usual Galilean 
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transformations of quantum mechanics. In ordinary quantum mechanics Schrodinger's equation has a time derivative 
that acts on the parameter in Uq . In the relational treatment the equation has a term involving pq instead of the time 
derivative. Notice that po is minus the total energy instead of just the "system energy" . Therefore the presence of the 
operator q° in Uq has the same effect in the relational treatment as the derivative with respect to the parameter has in 
ordinary quantum mechanics: they both induce a change in the energy of the system due to the boost, po ^ Po + (3t^- 
To study the changes in the conditional probability we go back in the derivation to equation (|51|l . 
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The value of the denominator actually does not change due to the boost, although its form changes. We will address 
this point later on. 

To understand the covariance it is convenient to commute Uq with the exponential; let us therefore analyze the 
product. 
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which can be shown using the Baker-Campbell-Hausdorff formula. 

We wish to commute Ug to the left. We start by noting that in the subspace ofci, tt, the variable q° behaves as an 
external parameter, as if it were a classical time t. Following the calculations of 8] for an ordinary quantum system 
one has, 
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which just states that the evolution of the Galilean transformed state should coincide with the Galilean transform of 
the original evolved state. That is, 
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We can therefore write for our system. 
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We still need to commute exp ( i (n — no) ) with Ug- Noting that the expression for Ug H65|I can be written 
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we see that only the first term in the exponential has a non-trivial commutator. 

To proceed we note that if one has two operators A, B such that [[A, -B], ^] — 0, one has that, 

If we now take A = a{po + £')^ and B ~ bcf we have the following identities, 
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Taking a — ^^"^""-^ and b = ^/3^ + /Stt and substituting po + E by po + E in (|68|l we finally have, 
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It is remarkable that all the terms involving i,„ax in (j68(l have cancelled with the terms stemming from the com- 
mutation we just did. We now can address the point we postponed before, namely the change in the denominator of 
the expression. Basically, a similar calculation to the one we just did starting from H45I) and performing the commu- 
tations shows that the denominator is actually invariant, using the cyclicity of the trace and the fact that unlike the 
numerator, it does not involve the projector on the (f> space. 

Now since 



UQp^^qoUG — P,j>-I3q0,q0, 

we therefore have, 

P' {(/> e A0|g" e Ai) = P (0 - /3g" e A0|g" e At) 
Which shows that the conditional probability is Galilean invariant. 



(76) 
(77) 



IV. CONCLUSIONS 



We have shown in a simple model that considering a quantum clock in quantum mechanics leads to a modification of 
Schrodinger equation, and that the resulting probabilities are Galilean invariant. Since the probabilities are invariant, 
then physical predictions from this framework will also be invariant. In particular the rate of decoherence predicted 
in 0, 1^ tlOj should be invariant. This is an interesting point since the rate of decoherence predicted is proportional 
to the difference of energies of states of the system in an basis of energy eigenstates. One could ask the question, 
how can this formula be Galilean invariant since the energy is not? The answer has to do with the fact that in order 
to have an energy basis as the one assumed in the calculation (with discrete spectrum) one has to consider systems 
analogous to a particle in a potential. In such systems, at least if they are isolated, the difference between energy 
levels is a Galilean invariant and therefore the decoherence rate is a Galilean invariant. 

It remains to be studied how the decoherence presented would transform under Lorentz transformations. Milburn 
|ll| recently studied decoherence in a Lorentz invariant setting and his treatment could provide a framework to analyze 
our proposal in some detail. This is more problematic, since we are considering corrections to quantum mechanics, 
and if one goes to the relativistic domain one first has to contend with the usual difficulties of defining a relativistic 
quantum mechanics. Although it appears that the use of a relational time could yield a well defined theory, the details 
remain to be studied. 
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